Magnetic field dependence of the superconducting gap node topology in 

non-centrosymmetric CePtsSi 
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The non-centrosymmetric superconductor CePtaSi is believed to have a line node in the energy 
gap arising from coexistence of s-wave and p-wave pairing. We show that a weak c-axis magnetic 
field will remove this line node, since it has no topological stability against time-reversal symmetry 
breaking perturbations. Conversely a field in the a — h plane is shown to remove the line node on 
some regions of the Fermi surface, while bifurcating the line node in other directions, resulting in 
two 'boomerang'-like shapes. These line node topological changes are predicted to be observable 
experimentally in the low temperature heat capacity. 



The existence of nodal points or lines in the energy gap 
of a superconductor or a Fermi superfliud, such as "^He, 
is one of the main characteristics of unconventional sym- 
metry pairing states. An important question is whether 
or not these nodes are accidental or are fuUy required 
by the fundamental symmetry of the pairing state. Sec- 
ondly one can ask whether or not the nodal points or lines 
are topologically stable, or whether they can in principle 
be destroyed by small perturbations. This question was 
discussed in a seminal paper by Volovikfl, in which he 
showed that point nodes, such as in ^He-A, are diabol- 
ical points characterized by a Berry phase and topolog- 
ical charge. Because of this property the nodal struc- 
ture cannot be removed by small perturbations, such as 
a magnetic field. Very recently Sato 2] and Volovik .3J 
have independently, shown that for line nodes of the gap 
function topological stability is only guaranteed if time 
reversal symmetry is preserved. The topological struc- 
ture of the nodal line can be classified as Z2 (i.e. the 
group, {0,1}, of integers modulo 2). Perturbations break- 
ing time reversal symmetry remove conservation of the 
relevant quantum numbers, and hence line nodes are not 
topologically stable against such perturbations. 

It is especially interesting to consider the gap nodal 
structure in the recently discovered non-centrosymmetric 
heavy fermion superconductor CcPt3Si4]. This mate- 
rial is found to have a spin-density wave Neel order 
at Tjv = 2.2K, corresponding to an effective moment 
O.16/iR/Ce0 and then superconductivity below about 
O.75K0. The superconductivity is expected to be nearly 
unique because this material is tetragonal with space 
group PAmm, having no center of inversion symmetry. 
The existence of such an inversion center is normally es- 
sential in the distinction between spin singlet and spin 
triplet superconductivity in systems where the spin-orbit 
coupling is large |^ 13, la, 0- Therefore, without an in- 
version center this system can be expected to possess 



simultaneous singlet and triplet pairing components on 
the sheets of its, non Kramers degenerate, Fermi surface. 
Thermal conductivity experiments suggest that the su- 
perconducting gap most probably has lines of nodes on 
the Fermi surface T^l, while NMR experiments [ill 
show evidence for crystal field excitations, as well as the 
formation of a heavy fermion metallic state. Surprisingly, 
1/TiT does not show the usual behaviour of other Ce 
based heavy fermion superconductors, but rather a weak 
Hebel-Slichter like peak below Tc. 

In this letter we consider the effect time reversal sym- 
metry breaking on the nodal topology in CePtaSi. Start- 
ing with the model gap nodal structure proposed by 
Hayashi et al.^^, we show that the line node can be 
completely removed by an arbitrarily weak c-axis mag- 
netic field. In contrast, for a — b plane oriented mag- 
netic fields we show that magnetic fields create additional 
lines of nodes, splitting the single degenerate node of the 
time-reversal symmetric ground state into two. More- 
over, we show that the double line of nodes do not ex- 



a)H = 



b) Hy^O 




node 




FIG. 1: (color online) Illustration of the influence of the ab- 
plane external magnetic field on the gap line nodes in CePtaSi. 
(a) zero-field line nodes for ±kz = const on one Rashba 
Fermi surface sheet; (b) Hy 7^ removes the node along 
the ±ky direction and bifurcates the node elsewhere, yield- 
ing "boomerang" -like nodal topology. 
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tend around the whole Fermi sphere but form a topo- 
logical "boomerang" -like structure, as shown in Fig. ^ 
These effects lead to a dramatic increase in the linear 
T coefficient in Cy /T, which should be visible experi- 
mentally. The appearance of these extra line nodes in 
a magnetic field and their large effect on the physical 
properties is surprising given the often stated argument 
that the non Kramers degenerate Fermi surface should 
suppress the effect of Zeeman interactions on the super- 
conducting state|l4|. However, a strong and unusual an- 
gle dependence to the paramagnetic susceptibility has 
already been predicted by Samokhin 1^ and Mineevfl^ 
has shown that the paramagnetic limit is expected to be 
strongly field- angle dependent. 

The superconducting state of CePtaSi has already been 
subjected to considerable experimental and theoretical 
effort. As well as the existence of the gap line node[lo|. 
a surprisingly high Hc2 has been observed of about 3.2T 
for H along (001) and 2.7T for H along (100) [13, exceed- 
ing the Pauli limiting field of Hp ^ 1.2T. In addition to 
the weak antiferromagnetism, with Q = (001/2) and an 
effective moment of about 0.16/xb per Ce atom, consider- 
ably less than the 2.54/iB/Ce expected for the Ce'^+ ion 
4/5/2 ground state, well defined crystal field excitations 
at ImeV and 24meV were also seen0. 

Theoretical models of the superconducting state in 
CePtaSi are based upon the existence of a Rashba type 
spin-orbit coupling 

Hso = " X! Sk-CTss'clsCks' (1) 
k.s,s' 
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as studied by Gor'kov and Rashba jlSj. Here gk = — g-k 
is a real pseudovector, by convention normalized so that 
idk) -F = 1 where (. . .)_f denotes an average over the Fermi 
surface, a = {(Tx, Uy, az) is the vector of Pauli matrices 
and c]^g is the usual Fermi field operator. This spin- 
orbit coupling removes the usual Kramers degeneracy 
between the two spin states at a given k, and leads to 
a quasi-particle dispersion Sk = Ck ± a|gk| splitting 
the Fermi surface into two sheets with a k dependent 
local spin orientation. The fuUy-relativistic Fermi sur- 
face calculated by Samokhin, Zijlstra and Bosejiflj shows 
three distinct sheets, each of which becomes split into two 
as a result of the spin-orbit interaction. The spin-orbit 
coupling breaks the parity symmetry P, and therefore 
mixes s ing let and triplet pairing states. A full symmetry 
analyis|2(J, |2l|, |23| shows that conventional s-wave pair- 
ing (Ak having full tetragonal symmetry) and a p-wave 
triplet pairing state of the form dk ~ i~ky, kx, 0) are able 
to coexist. If the microscopic pairing interaction strongly 
favors one type of pairing over the other then that one 
will be dominant, but in general both gap components 
are present on a single Fermi surface sheet. 



Following Hayashi et aL[l3j we take the pairing state 
as the combined s-wave and dk ~ {~ky, kx,0) state de- 
scribed above. Using the notation of Ref. the mean- 
field BCS Hamiltonian for this system in the presence of 
a magnetic Zeeman field can be written 
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where for brevity we write gk,2; + «gk,y = igwe^'^, with 
(j) = tan~^ (ky/kx). Similarly we write dk,a; -|- idk.j, = 
id\^e"^ and we have set Hx — 0. 

In zero external field the diagonalization of this Hamil- 
tonian yields the quasiparticle energy dispersion, 

i^i,2(k) - v/(ekTagk)2 + |dkTAk|2. (3) 

In the case where both singlet and triplet order parame- 
ters, Ak and dk, can be chosen as real, this is the spec- 
trum described by Hayashi et al^^. Assuming that the 
p-wave gap component is dominant so that Ak is smaller 
than the maximum value of dk on the Fermi surface, then 
one of the Rashba split sheets of the Fermi surface will 
have a line node, where Ak = dk and the other sheet will 
be nodeless. For a spherical Fermi surface and the sim- 



plest gap functions allowed by symmetry, Ak = A and 
dk = dQ(—ky,kx,0)/kp, the nodes are two horizonatal 
circles around the Fermi sphere at ±fcz = const, as illus- 
trated in Fig. n^a). As pointed out by Sato 2] the nodes 
have Z2 character. They can be removed by simply set- 
ting the relative phase between the two complex scalar 
order parameters Ak and dk non-zero, as one can readily 
see from Eq. |21 

Let us now consider what will happen to the lines of 
nodes when time-reversal symmetry is broken by apply- 
ing a Zeeman exchange field H. We assume that we 
are dealing with relatively weak field and therefore the 
Fulde-Ferrel-Larkin-Ovchinnikov statelH [H is not 
present. Let us first consider the case of a c-axis field. 
In this case H ± gk and the Hamiltonian (|2Jl can be 
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diagonalized to obtain 



and EsACk, Hz) = — i?i,2(k, i/^). One can readily see 
that this spectrum has no node, and so the gap node 
predicted by Hasyashi et aL is not topologicaUy stable 
against weak z-axis magnetic fields. 

To explain the lack of topological stability, consider 
the normal and superconducting state quasiparticle spec- 
trum near the two Rashba split sheets of the Fermi sur- 
face. In Fig[2Ia) where we show schematically the en- 
ergy dispersion Ek|| in a cut through the Fermi surface 
at kz — ±const. We have chosen the value of fc^ coincid- 
ing with the line node predicted by Hayashi et al., given 
by the zero of Eq.|3| On the left hand panel of Figl^Ja) we 
show the normal state electron (el, e2) and hole (hl,h2) 
quasiparticle dispersions near the two Rashba Fermi sur- 
face sheets. On the right hand side we show the cor- 
responding superconducting state energies from Eq. |31 
Clearly one sheet has a finite gap and the other has a 

a)H=0 




FIG. 2: (color online) Normal state (left) and superconduct- 
ing state (right) quasiparticle dispersions, and -Ek, for par- 
ticles and holes as a function of fe = (k'^ + ky)^^'^ for a fixed kz. 
From top to bottom, zero field H = 0, H J_ gk, and H |1 gk 
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node. Now consider the spectrum when a weak c-axis 
field is applied, as shown in Fig. [2Ib). In this case the 
spectrum becomes fully gapped, and the line node is re- 
moved. In effect the original line node is "accidental", 
arising from a particular cancellation between the s-wave 
and p-wave gap compoonents Ak and c?k at some partic- 
ular value of kz- When the symmetry breaking pertur- 
bation Hz is applied the crossing electron and hole-like 
levels h2 and e2 in Fig.[5fa), interact and so an avoided 
level crossing occurs, hence removing the line node. 

On the other hand consider the case of a a — 6 plane 
magnetic field. The two Fermi surface sheets will be 
perturbed differently depending upon whether H || gk 
or H _L gk. The latter case is essentially identical to 
Fig. Hfb), and so for these parts of the Fermi surface 
the gap node is removed. On the other hand, for the 
regions where H || gk the gap node is found to bifur- 
cate. The reason for this is illustrated in Fig.El^c). The 
field perturbs the two Rashba sheets oppositely, because 
they have spin components parallel to the applied field 
H. This can be seen in the normal state spectrum to the 
left side of Fig.Efc). The corresponding superconducting 
state spectrum is shown to the right in Fig. |2Ic). The 
four quasiparticle states no longer obey the symmetry 
E^ ^Ck, Hy) — —Ei^20^, Hy). A result is that one branch 
of the quasiparticle spectrum near the original line node 
now crosses zero twice, while the other remains non-zero. 
This implies that the gap node bifurcates. Note that the 
final spectrum can always be defined as positive, con- 
sistent with stability of the Fermi sea, after applying a 
suitable particle-hole transformation to the quasiparticle 
states. This sign change results in the bifurcated node, 
as shown in the right panel in Fig.lJfc). 

The final nodal topology for a — b oriented magnetic 
fields is shown in Fig. ^h). For k vectors with x — y 
components parallel to H the gap node is removed, while 
for k vectors with x — y components perpendicular to H 
the gap node bifurcates. Overall this leads to the unusual 
"bomerang" shaped gap line nodes as shown in Fig.^b). 
This topological structure is consistent with the general 
Z2 topologigal arguments of Sato0 and VolovikQ. The 
Z2 gap node has two topological invariants, which in the 
presence of time reversal symmetry are independently 
conserved, leading to stability of the line node. But the 
symmetry breaking perturbation H removes the stabil- 
ity, leading to either node removal or node bifurcation 
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FIG. 3: (color online) Calculated temperature dependence 
of the specific heat with and without Zeeman field along z- 
direction(a) and along j/-direction(b). The inset shows the 
low-temperature region for H-^ — and = 3T. Here, we 
assume the BCS-like temperature dependence of the both p- 
wave and s-wave components of the superconducting gap with 
dk(0) = 0.2meV and Ak(0) = 0.25Ap(0) which gives about 
the BCS ratio between the magnitude of the superconduct- 
ing gap and Tc = 0.75K. Following experimental findings we 
adopt Ef ~ 0.012eV and estimate kp ~ 0.5 which are 
typical values for the heavy-fermion metals. We further ap- 
proximate the Rashba splitting of the bands a ~ 1 x 10~^ eV 

A. 

depending on the field orientation. 

Experimentally this nodal change should be observable 
in unusual field dependence of various low temperature 
properties. In CePtaSi NMR experiments show both a 
Hebel-Slichter peak, a signature of the nodeless s-wave 
superconductor, and power-law behavior of the spin- 
lattice relaxation at low temperatures consistent with 
nodal structure of the superconducting gap [H [13. The 
predicted gap topology in Fig. ^ would imply a change in 
the low temperature power law to exponential for c-axis 
fields, and a changed prefactor of the low temperature 
power law for a — 6 oriented fields. Similarly the changes 
in the line node topology upon applied magnetic field re- 
sults in the change of the thermodynamic characteristics 
such as heat capacity, Cv/T. In FigOwe show our re- 
sults for the temperature dependence of the specific heat, 
Cv/T for various Zeeman fields. For c-axis fields the low 
temperature behavior changes from linear to exponential, 
but the actual total change is small for physically relevant 
values of , as shown in Fig|3fa) . In contrast for a — b 
plane fields the formation of the "boomerang" -like struc- 
ture in the line nodes nodes both reduces the jump in the 
specific heat at Tc, and leads to a noticeable increase in 
Cv /T at low temperatures. Furthermore, for larger fields 
we also find that an even more substantial reduction of 
the jump occurs. This happens because for large enough 
exchange fields a new pair of gap nodes develops also on 
the second Rashba Fermi surface sheet. This is because 
for larger fields the lower quasi-particle state E2 (k) , also 
crosses zero as is evident from Figl^Jc). Most interest- 
ingly we find that these effects develop at small enough 
fields in CePtaSi to be observed experimentally. 



In conclusion, we have shown that Zeeman field in- 
duces an interesting and novel phenomena in the the non- 
centrosymmetric heavy fermion superconductor CePtaSi. 
A c-axis field is shown to remove the "accidental" line 
node, and change the specific heat, Cy/T from linear 
to exponential at low temperatures. A field in the a — b 
plane is shown to spit the nodes into two 'boomerang'-like 
structures which will reduce the jump in the Cy/T and 
enhance the value of Cv/T at zero temperature. These 
surprising changes result from the changes of the Z2 topo- 
logical quantum numbers by the time reversal symmetry 
breaking perturbation, consistent with the arguments of 
Sato0 and Volovik 0. 
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